CONVERGENCE OF EQUILIBRIA OF THIN ELASTIC PLATES 
UNDER PHYSICAL GROWTH CONDITIONS FOR THE 
ENERGY DENSITY 



MARIA GIOVANNA MORA AND LUCIA SCARDIA 

Abstract. The asymptotic behaviour of the equilibrium configurations of a 
thin clastic plate is studied, as the thickness h of the plate goes to zero. More 
precisely, it is shown that critical points of the nonlinear elastic functional E h , 
whose energies (per unit thickness) are bounded by Ch 4 , converge to critical 
points of the T-limit of h~ 4 £ h . This is proved under the physical assumption 
that the energy density W(F) blows up as det F — > 0. 



1. Introduction 

A thin plate is a three-dimensional body, occupying in a reference configuration a 
region of the form f2/j := 5x(— ^, ^), where the mid-surface S is a bounded domain 
of M 2 and the small parameter h > measures the thickness of the plate. 

The elastic behaviour of such bodies is classically described by means of two- 
dimensional models, which are easier to handle both from an analytical and a 
computational viewpoint than their three-dimensional counterparts. There exists 
a large variety of such theories in the literature (see [5, 12] for a survey). However, 
as their derivation is usually based on a priori assumptions on the form of relevant 
deformations, their rigorous range of validity is typically not clear. A fundamen- 
tal question in elasticity is thus to justify rigorously lower dimensional models in 
relation to the three-dimensional theory. 

Recently, a novel variational approach through T-convergence has led to the 
rigorous derivation of a hierarchy of limiting theories. Among other features, it 
ensures the convergence of three-dimensional minimizers to minimizers of suitable 
lower dimensional limit energies. 

In this paper we discuss the convergence of (possibly non-minimizing) stationary 
points of the three-dimensional elastic energy, assuming physical growth conditions 
for the energy density. 

We first review the main results of the variational approach. Given a thin plate 
Qh, the starting point of the variational analysis is the three-dimensional nonlinear 
elastic energy (scaled by unit thickness) £ h (w,$lh) associated to a deformation w 
of the plate. The limiting behaviour of £ h as the thickness of the plate tends to 
zero, can be described by the T-limit Ip of the functionals 

h-P£ h (;Q h ), 
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as h — ► 0, for a given scaling (3 > 0. As mentioned above, this implies, roughly 
speaking, convergence of minimizers w h of £ h (-,fth) (subject to applied forces or 
boundary conditions) to minimizers of the two-dimensional energy Tp, provided 
£ h (w h , f2/j) < ChP . For the definition and main properties of T-convergence we 
refer to the monographs [4, 7]. 

In this setting, T-convergence was first proved by Le Dret and Raoult in [10] for 
the scaling (3 = 0. This led to a rigorous justification of the nonlinear membrane 
theory. In the seminal papers [8, 9] Friesecke, James, and Muller established T- 
convergence for all (3 > 2. The scaling (3 = 2 corresponds in the limit to the 
Kirchhoff plate theory, while j3 = 4 to the von Kdrmdn plate theory. For (3 > 4 the 
usual linear theory is derived, while the intermediate scalings 2 < (3 < 4 relate to 
a linear theory with constraints. The case of < (3 < 5/3 was recently solved by 
Conti and Maggi [6]. The regime 5/3 < (3 < 2 remains open and is conjectured to 
be relevant for crumpling of elastic sheets. Analogous results have been proved for 
thin rods in [1, 15, 16]. 

The intent of this paper is to investigate the convergence of stationary points 
of the three-dimensional nonlinear elastic energy (subject to applied forces and 
boundary conditions) to stationary points of the r-limit functional. The first result 
concerning convergence of equilibria for thin bodies has been shown in [17], in the 
case of a thin strip and for the scaling (3 = 2. This work has been extended in [18] 
to the case of a thin rod in the regime (3 = 2, and then in [19] to a thin plate in the 
von Karman regime (3 = 4 (see also [11] for an extension to thin shells). A crucial 
assumption in all these papers is that the elastic energy density W is differentiablc 
everywhere and its derivative satisfies a linear growth condition. This bound on 
DW is unsatisfactory, as it prevents the blow-up of W(F) as the determinant of the 
deformation gradient F tends to zero (which corresponds to a strong compression 
of the body). We point out that, instead, the results in [8, 9], as well as the ones in 
[15, 16], do not require any bound from above on W . On the other hand, without 
assuming a linear growth condition on DW, it is not even clear to which extent 
minimizers satisfy the Euler-Lagrange equations in the conventional form (see (2.10) 
below) . 

A growth condition on W, which is compatible with the blow-up condition as 
dct F -> is: 

\DW{F)F T \ < k(W(F) + 1) (1.1) 

for every F with dctF > 0. In [2, 3] Ball has shown that, under assumption (1.1), 
it is possible to derive an alternative first-order necessary condition for minimizers 
(Theorem 2.1). When minimizers are invertible, this condition corresponds to the 
equilibrium equation for the Cauchy stress tensor (Remark 2.2). 

In this paper we focus on the scaling (3 = 4 and we extend the analysis of [19] 
to the case of an elastic energy density W satisfying the physical growth condi- 
tion (1.1). More precisely, we call a deformation a stationary point of the three- 
dimensional energy if it satisfies the first-order necessary condition introduced by 
Ball in [2, 3] (Definition 2.3). In Theorem 3.1 we prove that any sequence of sta- 
tionary points w of the three-dimensional energy, satisfying £ h (w h , < Ch A , 
converges to a stationary point of the von Karman functional (i.e., to a solution 
of the classical Euler-Lagrange equations of the limiting functional). This is the 
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first result of convergence of equilibria for thin plates compatible with the physical 
requirement that W(F) — > +00 as dctF — > 0. 

A key ingredient in the proof of our main result is the quantitative rigidity 
estimate proved by Friesecke, James, and Miiller in [8, Theorem 3.1]. It is first 
used to deduce compactness of sequences of stationary points from the bound on 
the elastic part of the energy, and then to define suitable strain-like and stress-like 
variables G h and E h (see (4.9) and (4.13)). While in [19] the 1? bound on the 
strains G h (which follows directly from the rigidity estimate) implies immediately 
an analogous bound for the stresses E h , the present case exhibits an additional 
difficulty. Indeed, in our setting the stresses E h turn out to be naturally defined as 

E h = ^DWild + h 2 G h )(Id + h 2 G h ) T (1.2) 
h 

(E h can be interpreted as a sort of Cauchy stress tensor, read in the undeformed 
configuration, see also Remark 2.2). Hence, using the growth condition (1.1) and 
the bound on the elastic energy we can only deduce weak compactness of E h in L 1 
and this convergence is not enough to pass to the limit in the three-dimensional 
Euler-Lagrange equations (see Steps 2 - 3 of the proof and the discussion therein). 

This difficulty is overcome by identifying a sequence of measurable sets Bh , which 
converge in measure to the whole set fi := Sx(—h,^) and satisfy the following 
properties. On the remainder in the first order Taylor expansion of DW around 
the identity is uniformly controlled with respect to h, so that one can deduce an L 2 
bound for E h from (1.2) and from the L 2 bound on G h . On the complement of Bh 
one can use the growth condition (1.1) to show that the contribution of E h on this 
set is negligible at the limit in the L 1 norm. This mixed type of convergence of the 
stresses is then shown to be sufficient to pass to the limit in the three-dimensional 
Euler-Lagrange equations. 

Our main result also applies to the regime (3 > 4. More precisely, in Theorem 3.1 
we also show that any sequence of stationary points w h of the three-dimensional 
energy, satisfying £ h (w h ,Vl h ) < Ch 13 with (3 > 4, converges to a stationary point of 
the functional of the linear plate theory. 

Convergence results for thin plates in the Kirchhoff regime (3 = 2 and in the in- 
termediate scalings 2 < f3 < 4 are still open, even under the simplifying assumption 
of linear growth of DW. The additional difficulties in the analysis of these regimes 
are due to the weaker compactness properties arising from the rigidity theorem 
and to the presence, in the limiting model, of a nonlinear or geometrically linear 
isometry constraint. 

The plan of the paper is the following. In Section 2 we describe the setting of 
the problem and we discuss the first order necessary condition by Ball. Section 3 
contains the statement of the main result, which is proved in Section 4. 

2. Setting of the problem 

We consider a thin plate, whose reference configuration is given by the set Clh = 
Sx ( — -|, -|), where S c K 2 is a bounded domain with Lipschitz boundary and h > 0. 

Deformations of the plate are described by maps w : flh — > K 3 , which are assumed 
to belong to the space ff 1 (0? l ;R 3 ). Moreover, we require the deformations w to 
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satisfy the boundary condition 

w(z)=z for every z G Tx(-|, §), (2.1) 

where T is a (non-empty) relatively open subset of dS. 

To any deformation w G ^(Sl^K 3 ) we associate the total energy (per unit 
thickness) defined as 

F h (w) = \f W(Vw)dz-y f f h -wdz, (2.2) 
h Jn h h J nh 

where f h G L 2 (f2/,.; R 3 ) is the density of a body force applied to f^. 

On the stored-energy density W: M 3x3 — > [0, +00] we require the following 
asssumptions: 

W is of class C 1 on M 3x3 ; (2.3) 
VK(F) = +00 ifdetF<0, -> +00 asdetF->-0 + ; (2.4) 

W(iLF) = W(.F) for every R G 5*0(3), F G M 3x3 (frame indifference). (2.5) 

Here M 3x3 denotes the set of matrices F G M 3x3 with dctF > 0, while SO (3) 
denotes the set of proper rotations {R G M 3x3 : R T R = Id, det R = 1}. Condition 
(2.4) is related to the physical requirements of non-interpenetration of matter and 
preservation of orientation. It ensures local invertibility of C 1 deformations with 
finite energy. 

We also require W to have a single well at 50(3), namely 

W = on 50(3); (2.6) 
W{F) > Odist 2 (F, 50(3)); (2.7) 
W is of class C 2 in a ^-neighbourhood of 50(3). (2.8) 
Finally, we assume the following growth condition: 

\DW{F)F T \ < k(W(F) + 1) for every F G M 3x3 . (2.9) 

This is a mild growth condition on W , introduced by Ball in [2, 3], which is com- 
patible with the physical requirement (2.4), but is nevertheless sufficient to derive 
a first-order condition for minimizcrs of T h . In fact, by performing external varia- 
tions w + £<j) of a minimizer w, one is formally led to the Euler-Lagrange equations 
in the conventional form 

/ DW{Vw) -V(j)dz = f h -(j)dz V^> smooth with ^>| rx( _h ) h ) = 0. (2.10) 

To justify rigorously this derivation, one has to require that either DW is Lipschitz 
continuous or the minimizer w belongs to W 1,ao and satisfies a stronger orientation 
preserving condition, namely det Ww > c > a.c. in Cl^. However, none of these 
assumptions is satisfactory: the Lipschitz continuity of DW is incompatible with 
(2.4), while there may exist minimizcrs that do not belong to W 1 ' 00 or do not satisfy 
the stronger orientation preserving condition (see the discussion in [3, Section 2.4]). 

If instead condition (2.9) is assumed, then it is possible to derive an alternative 
equilibrium equation for minimizers. More precisely, by considering variations of 
the form w + £ 4> o w one can deduce the following condition. 
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Theorem 2.1 ([3, Theorem 2.4]). Assume thatW satisfies (2.3), (2.4), and (2.9). 
Let U C M 3 be a bounded domain with a Lipschitz boundary dU = dU\ U dUi U 
N, where dU\, dU2 are disjoint relatively open subsets of dU and N has zero 
two-dimensional measure. Let w G H 1 / 2 {dU; R 3 ) and f G L 2 ([/;R 3 ). Let w G 
ii^f/jR 3 ) 6e a /ocaZ minimizer of the functional 

T(w) := [ W{Vw)dz- [ f-wdz 
Ju Ju 

subject to the boundary condition w — w on dU\, that is, there exists e > such 
that J-(w) < J~(v) for every v G H l {U\ R ) satisfying \\v — < £ and v = w on 

dU x . Then 

I DW{\7w)(\7w) T :V<f){w)dz = f f-<f>{w)dz (2.11) 
Ju Ju 
for all <j) G Cl (M 3 ; M 3 ) such that <fi ° w — on dU\ in the sense of trace. 

In the theorem above and in what follows, given a subset U of l n we denote 
by C£(U) the space of functions of class C k that are bounded in U, with bounded 
derivatives up to the fe-th order. We also stress that in (2.11) the term V</>(w) 
denotes the gradient of <j) computed at the point w{z). 

Remark 2.2. Under the assumptions of Theorem 2.1, if in addition w is a smooth 
homcomorphism of U onto U' := w(U), then equation (2.11) reduces by means of 
a change of variables to 

f T(w~ 1 (x)):V(j>{x)dx = [ } '(uT 1 (x)) ■ <j>(x) dx 

Jw(U) Jw{U) 

for all (f> G C 1 (M 3 ; R 3 ) such that <p\ w (dUi) — 0- I n the formula above T is the Cauchy 
stress tensor: 

T(z) = (det Vw{z))- 1 DW{Vw{z)){y W {z)) T , zeU 

and / = (det Vw)~ 1 f (see [3, Theorem 2.6]). In other words, Theorem 2.1 asserts 
that the equilibrium equations are satisfied in the deformed configuration. 

In our setting it is natural to assume (2.11) as definition of stationary points 
of T h . Our aim is to analyse their limit behaviour, as the thickness h goes to 0. 
To do so, it is convenient to perform a change of variables and to reduce to a 
fixed domain independent of h. Thus, we consider the scaling (z',z 3 ) = (x',hx^), 
Vh = (V, 7^3), y(x) = w(z), and g h (x) = f h (z), and we introduce the functional 

J h (y)=J*(w)= [ W(V h y)dx- [ g h -ydx, (2.12) 
Jn Jn 

where il = Sx(—^,h) and the scaled deformation y G if 1 (il;R 3 ) satisfies the 
boundary condition 

y(x) = (a/, hx 3 ) for every x = (x', X3) G Fx(— |, |). (2.13) 

According to Theorem 2.1, we give the following definition. 

Definition 2.3. We say that a deformation y G i? 1 (J7;R 3 ) is a stationary point of 
J h , subject to clamped boundary conditions on Tx(— ^, |), if y(x) = (x', hxj) for 
every x G Tx( — i, |) and the following equation is satisfied: 

/ DW(V h y){V h y) T :V^{y)dx= [ g h ■ 4>(y) dx 
Jn Jn 
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for all (j> G Cl(M. 3 ;R 3 ) satisfying <f>(x',hx 3 ) = for every x G rx(-|, \). 

Remark 2.4. In [3] Ball has shown that, if W satisfies the growth condition 

\F T DW(F)\ < k(W(F) + 1) for every F G M 3x3 

(which implies, but is not equivalent to (2.9), sec [3, Proposition 2.3]), then local 
minimizcrs of T satisfy the equation 

(W(Vw)Id- (Vw) T DW(\7w)) : Vcpdz = [ (\7w) T f ■ <j)dz 



for all 4> G Cq (U; R 3 ). This equation is obtained by performing internal variations of 
the form woip e: with ip~ 1 (x) — x + s<fi(x), and can be viewed as a multi-dimensional 
version of the classical Du Bois- Raymond equation of the one-dimensional calculus 
of variations. To the purpose of our analysis the use of this equilibrium equation in 
place of (2.11) seems to be less convenient. Indeed, the requirement of zero bound- 
ary values for the test functions suggests that the equation does not provide precise 
information about the boundary behaviour of the limiting quantities. Moreover, it 
imposes a severe restriction on the choice of admissible test functions. 

Remark 2.5. If W satisfies (2.9), then W has polynomial growth, that is, there 
exists s > such that 

W(F) < C{\F\ S + IF^I 8 ) for all F G M^ x3 

(see [3, Proposition 2.7]). In particular, examples of functions satisfying (2.3)-(2.9) 
are: 

W{F) = \{F T F) 1/2 -Id| 2 + |logdetF| p for F G M^ 3 , 

or 

W(F) = \(F T F)V 2 -Id] 2 -, 
where p > 1 and W is intended to be +oo if det F < 



1 -1 



for F G M 3x3 , 



3. Statement of the main result 

In this paper we focus on the asymptotic study of stationary points y h of J h 
(according to Definition 2.3) with elastic energy (per unit thickness) of order h 13 
with ft > 4, that is, 

/ W(V h y h )dx < Ch , /3 > 4. (3.1) 

For simplicity we assume that the body forces g h are independent of the variable 
£3 and normal to the mid-surface of the plate; more precisely, we assume g (x) = 
h^ +2 ^ 2 g(x')es, where g G L 2 (S) is given. The scaling /i(^+ 2 )/ 2 of the normal force 
ensures consistency with the elastic energy scaling (3.1). 

In [9] Friesecke, James, and Midler have identified the limit of the functionals 
j , in the sense of T-convergence, under the assumptions (2.5)-(2.8). For [3 = 4 
the P-limit ^7 v k can be expressed in terms of the averaged in-plane and out-of-plane 
displacements u and v (see (3.13)) and is given by 



J v k(u,v) = l vK {u,v) ~ 
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where, for u £ H 1 {S] R 2 ) and v e H 2 (S), the von Karman functional 2" v k is defined 
as 

I vK (u,v) = ^J Q2(symV'u + ^V'v®Vv)dx' + — J Q 2 ((V) 2 v) dx' . (3.3) 

Here Q 2 is a quadratic form that can be computed from the linearization of W at the 
identity. More precisely, we consider the quadratic form Q$(F) = D 2 W(ld)F : F on 
M 3x3 and define the quadratic form Q 2 on M 2x2 through the following minimization 
procedure: 

Q 2 (G) = C 2 G:G := min Q 3 (F), (3.4) 

F" — G 

where F" denotes the 2x2 submatrix given by F/j = Fij, 1 < i,j < 2. 

For (3 > 4 the T-limit J\[ n depends only on the averaged out-of-planc displace- 
ment v and is given by 



J\in(v) = I\in{v ) - / gvdx', (3.5) 
Js 

where Ty m is the functional of linear plate theory, defined as 

Ilin(v) = Ya f Q2{{V'Yv)dx'. (3.6) 

J s 

for every v G H 2 (S). 

The T-convergence result guarantees, in particular, that given a minimizing se- 
quence y h satisfying 

]imsup^(j h (y h )-MJ h ) =0, 

the averaged in-plane and out-of-plane displacements associated with y h converge 
to a minimizer (u, v) of J7vK if /3 = 4. If (3 > 4, they converge to a pair of the form 
(0,v), where v is a minimizer of J\ m . 

To set the stage for our result on the convergence of equilibria, we derive the 
Euler-Lagrange equations for a minimizer (it, v) of J v k- First of all, from the 
clamped boundary conditions (2.13) it follows that the limiting displacement (u, v) 
satisfies 

u{x') = and v(x') = 0, W'v(x') = for every x' E T. (3.7) 

By computing the respective variations of J7vK in u and v we obtain the following 
Euler-Lagrange equations in weak form: 

[ (^(symV'u + ]-Vv ® V'v) : {V'v ® VV) 

Js V 2 (3.8) 

+ ^C 2 ((V) 2 v):(V) 2 p-g^j dx' = 

for every ip e H 2 (S) with ^| r = 0, VVIr = 0, and 

[ C 2 (sym V'u + ]-Vv® V'v) : V'lp dx' = (3.9) 
Js 2 

for every tp e iJ x (5;R 2 ) with V'|r = 0. 
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In the case of the linear functional J7n n the limit displacement v satisfies the 
boundary conditions 

v(x') = 0, V'v(x') = for every x' G T (3.10) 

and the Euler-Lagrange equations are given by 

J (^C 2 ((V') 2 v) : (V')V - 9?) dx' = (3.11) 

for every <p G H 2 (S) with <p\ T = 0, VVIr = 0, 

From now on we will adopt the notation y = (y',y 3 ). 
The main result of the paper is the following. 

Theorem 3.1. Assume that the energy density W satisfies (2.3) -(2.9). Let j3 > 4. 

Let (y h ) be a sequence of stationary points of J h according to Definition 2.3, with 
g h (x) = h^ +2 '' 2 g(x')e 3 . Assume further that 

I W(V h y h )dx < Ch' 3 . (3.12) 
Jq 

Set 

1 

u\x>) := f 2 ((y h y(x',x 3 )-x')dx 3 , 

h x (3.13) 
h (P-2)/2 / t y 3 (x',x 3 )dx 3 

for every x' G S. Then the following assertions hold. 

(i) (von Kdrmdn regime) Assume j3 = 4. Then, there exist u G -ff 1 (S';K 2 ) and 
v G H 2 (S) such that, up to subsequences, 

u h ^u weakly m i7 1 (5;R 2 ) (3.14) 

and 

v h w strongly m H 1 (S), (3.15) 

as ft, — > 0, and £/ie Zimit displacement (u,v) solves (3.8)- (3.9), and satisfies 
the boundary conditions (3.7). 

(ii) (linear regime) Assume (3 > 4. Then, (3.14) and (3.15) hold with u = 0, and 
the limit displacement v solves (3.11) and satisfies the boundary conditions 
(3.10). 

Remark 3.2. If y h is a sequence of minimizers of J h with g h (x) = h^ +2 ^ 2 g{x')e 3 , 
then condition (3.12) is automatically satisfied. This can be proved by means of 
a Poincare-like inequality related to the rigidity theorem by Friesecke, James, and 
Miillcr (see the proof of [9, Theorem 2, part hi]). 

Remark 3.3. In [13] Mielkc used a centre manifold approach to compare solutions 
in a thin strip to a one-dimensional problem. This method works already for finite 
h, but it requires that the nonlinear strain (V; l j/) T V? l w is close to the identity in 
L°°. Applied forces g are also difficult to include. We also mention a more recent 
result by Monncau [14], based on a careful use of the implicit function theorem. 
He shows that, starting from a sufficiently smooth solution of the von Karman 
equations, there exists a nearby solution of the three-dimensional problem. 
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4. Proof of Theorem 3.1 
This section is devoted to the proof of the main result of the paper. 
Proof of Theorem 3.1. Let (3 > 4. For notational convenience we set 

a:= G9 + 2)/2, 

so that a > 3. Let (y h ) be a sequence of stationary points of J h , i.e., suppose that 

DW(V h y h )(V h y h f :V<f>(y h )dx = f h a ge 3 ■ <f>(y h ) dx (4.1) 
a Jn 

for all 4> G Cl(R 3 ;R 3 ) satisfying (f)(x',hx 3 ) = for every x € rx(-§, |). Further- 
more, assume that condition (3.12) is fulfilled. 

Step 1. Decomposition of the deformation gradients in rotation and strain. The 
energy bound (3.12) and the coercivity condition (2.7) imply that 



/ dist 2 (V h y h ,SO (3)) dx < Ch 
Jn 



2a-2 



Owing to the rigidity estimate [8, Theorem 3.1], this bound guarantees the existence 
of a sequence of smooth rotations R h , whose I? distance from V? l y' 1 is of order h a ~ 1 . 
A careful analysis of the increment of R h in neighbouring squares of side h shows 
that the gradient of R h is well controlled in terms of h. From this it follows that 
V^y' 1 must converge to a constant rotation (namely, the identity, because of the 
boundary condition) and that the in-plane and out-of-plane displacements satisfy 
the compactness properties (3.14) and (3.15), respectively. 

More precisely, arguing as in [9, Theorem 6 and Lemma 1], one can construct a 
sequence (R h ) C C°°(S; M 3x3 ) such that R h {x') G SO{3) for every x' G S and 

||V fc y h -iJ h || ia <C7i 0, - 1 I (4.2) 

l|V'i?' l || L 2 < C7i Q - 2 , (4.3) 

\\R h -Id|| L 2 < Ch a ~ 2 . (4.4) 

From (4.2) and (4.4) it follows that X7 h y h converge to Id strongly in L 2 (f2; M 3x3 ); 
in particular, Vy h — > diag{l, 1, 0} strongly in L 2 (f2;M 3x3 ). Therefore, by the 
boundary condition y h (x' , X3) = (x' , hx 3 ) for every x G Lx(— i, i) and the Poincare 
inequality, we have that 

y h -> (x', 0) strongly in H 1 ^; R 3 ). (4.5) 

By [9, Lemma 1] there exist u £ ^{S^R 2 ) and v G H 2 (S) such that (3.14) and 
(3.15) hold true, up to subsequences. From the boundary condition satisfied by y h 
we obtain immediately that u(x') = and v(x') = for every x' G T. Moreover, 
by [9, Corollary 1] the first moment £ /l of the in-plane displacement satisfies 

£V):=77^T r^y^V,^)-^')^ - -T^' v weakly in H 1 (S;R 2 ). 

As £ h = on T for every h, this implies V'w = on T. Finally, [9, Lemma 1] 
guarantees the following convergence properties for R h : 

r>h tj 

A h :=-f^- 2 >■ A:=-(V'v,0)(x)e 3 + e 3 ®(V'v,0) in H 1 (S; M 3x3 ) (4.6) 
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and 

s y m -^r^l" ini 9 (5;M 3x3 ), Vg < oo. (4.7) 

In particular, by the Poincare-Wirtingcr inequality and the equations (4.2) and 
(4.7), we obtain 

fl- x ^ h ^\h\\ L2 <c||^£ _i|| l2 <Ch a ~\ (4.8) 

The bound (4.2) suggests the following decomposition for the deformation gra- 
dients: 

V h y h = R h (ld + h 01 - 1 ^ 1 ). (4.9) 

By (4.2) the G h : Q -> M 3x3 are bounded in L 2 (VL; M 3x3 ). Thus, up to subse- 
quences, G h G weakly in L 2 (ft;M 3x3 ) for some G 6 L 2 (f7;M 3x3 ). By [9, 
Lemma 2] the limiting strain G satisfies 

G"(x', x 3 ) = G (x') ~ x 3 (V) 2 v, (4.10) 

where 

symG = symV'it + ^V'u ® Vv if a = 3, (4-11) 

symG = symV'u if a > 3. (4-12) 

We recall that G" denotes the 2x2 submatrix G"j = Gij, 1 < i,j < 2. 

S'tep ,2. Estimate for the scaled stress. Let E'' 1 : — > M 3x3 be the scaled stress 
defined by 

E h := — l_DW(Id + /t° _1 G fc )(Id + h a - 1 G h ) T . (4.13) 

Notice that -E^ is symmetric, due to the frame indifference of W. Moreover, the 
following estimate holds true: 

Indeed, if /i a " 1 |G ,l | < 5/2, where 5 is the size of the neighbourhood in (2.8), then 

DW (Id + h^G' 1 ) = h a - 1 D 2 W{F h )G h , 

for some matrix F h G M 3x3 with - Id| < 5/2. As D 2 W is bounded in this set, 
we deduce that 

\DW "(Id + h a ~ 1 G h )\ < Ch a - 1 \G h \, 

which implies 

\E h \ <C\G h \ + Ch a - 1 \G h \ 2 <C(l + 5)\G h \. 
If h a - 1 \G h \ > 5/2, by (2.9) we have 

1^1 < j^k(W(Id + h^G h ) + 1) < k W(ld + h^) + 2fc |Gfc| _ 

We notice that we are allowed to use the bound (2.9), as W(\7 h y h ) is finite a.e. 
in Q by (3.12), hence det V h y h = dct(Id + h a ~ x G h ) > a.e. in fi. This concludes 
the proof of (4.14). 
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Step 3. Convergenge properties of the scaled stress. By the decomposition (4.9) 
and the frame indifference of W , we obtain 

DW(V h y h )(V h y h f = R h DW{\& + h a - l G h ){\& + hf*' 1 G h ) T {R h ) T 
= h a - 1 R h E h {R h ) T . 

Thus, in terms of the stresses E h the Euler-Lagrange equations (4.1) can be writ- 
ten as 

f R h E h (R h ) T ' :V<t>{y h )dx= f hge 3 ■ <j>{y h ) dx (4.15) 
Jn Jn 

for all (f> G Cl(R 3 ;R 3 ) satisfying hx 3 ) = for every x G rx(-±, \). 

In order to pass to the limit in (4.15) we are interested in studying the conver- 
gence properties of the scaled stresses E h . 

By (4.14), (3.12) and the fact that the G h arc bounded in L 2 (fl; M 3x3 ), we 
deduce that for every measurable set A C £1 

\E»\dx < Gf W ^ h T Gh) dx + cf\G^dx 

J A J A 

< Ch *- 1 + C\A\ 1/2 . (4.16) 

This bound ensures that the scaled stresses E h are bounded and equi-integrable in 
L 1 (fi;M 3x3 ). Therefore, by the Dunford-Pcttis theorem 

E h ->> E weakly in L^Q; M 3x3 ) (4.17) 

for some E G L 1 (ri;M 3x3 ). In particular, since E h is symmetric, also E is sym- 
metric. 

One can immediately realize that weak convergence of E h in L 1 is not enough 
to pass to the limit in (4.15). This is due to the fact that, for instance, one cannot 
guarantee uniform convergence of the term V(j>(y h ) (recall that for y h we have the 
convergence (4.5)). Therefore, some more refined convergence properties for E h 
are needed. In particular, we shall identify a sequence of sets Bh, whose measures 
converge to the measure of fi (and therefore, on fl \ Bh the sequence E h converges 
to strongly in L 1 by (4.16)), and such that on Bh the sequence E h is weakly 
compact in L 2 . Using the C\ regularity of test functions, we shall show that this 
mixed type of convergence is sufficient to derive the limit equations. 

Let B h := {x G Q : h a - 1 -^\G h {x)\ < 1}, with 7 G (0, a — 2), and let X h be its 
characteristic function. Notice that 

|n\fl fc |< / h a - 1 -~<\G h \dx<Ch a - 1 -~ / \Sl\B h \ 1 / 2 \\G h \\ L 2, 
Jn\B h 

hence 

\n\B h \ < Ch 2 ^- 1 -^. (4.18) 

This implies in particular that \h converges to 1 in measure and thus, XhG h con- 
verges to G weakly in L 2 (f2;M 3x3 ). 
From (4.16) and (4.18) it follows that 

\E h \dx < Ch 01 - 1 ^, (4.19) 

Q\B h 

hence 

(l-Xh)E h ^0 strongly in L 1 (fi;M 3x3 ). (4.20) 



12 



M.G. MORA AND L. SCARDIA 



On the set Bh we have a uniform control of the term h a ~ 1 G h , so that we can 
deduce weak convergence of XhE h in L 2 (f2; M 3x3 ) from the weak convergence of 
G h simply by Taylor expansion. More precisely, let C be the linear operator defined 
by C := D 2 W(ld). We claim that 

XhE h CG weakly in L 2 {fl; M 3x3 ). (4.21) 

We note that, as R h converges boundedly in measure to Id, the claim implies 
also that XhR h E h converges to CG weakly in L 2 (fl; M 3x3 ). This remark will be 
repeatedly used in the next steps of the proof. 
By Taylor expansion we have 

£>W(Id + h^G 11 ) = h a - x CG h + r ] {h a ~ 1 G h ) 1 

where the remainder rj satisfies r/(F)/\F\ — > 0, as \F\ — > 0. This identity leads to 
the following decomposition of XhE h : 

XhE h ^ X h 1 ^{h a - 1 CG h +f ] {h a - 1 G h )){lA + h a - 1 G h ) T 



= Xh CG h + Xhh a ~ 1 £G h (G h ) T + Xh ^l^ + XhV{h"- 1 G h )(G h ) T . (4.22) 

To prove the claim (4.21) we analyse carefully each term on the right-hand side of 
(4.22). The weak convergence of XhG h to G in L 2 {Q; M 3x3 ) and the linearity of C 
yield 

£(XhG h ) CG weakly in L 2 (ft;M 3x3 ). (4.23) 
The second term in the right-hand side of (4.22) can be estimated as follows: 

\ X hh a - 1 CG h (G h f\ < XhCh^lG' 1 ] 2 < Ch^\G h \. 

Therefore, it converges to zero strongly in L 2 (fl; M 3x3 ) by the L 2 bound of the G h . 
As for the third term in (4.22), we have the following bound: 

77(/i Q - 1 G' 1 ) 



where for every t > we have set 



< oj(h^) \G h \ 



u>(t) :=sup|-j-^- :\A\<t 

Since w(t) — > for t — > + , we can conclude as before that Xh r l{h a ~ 1 G h ) //i™ -1 
converges to zero strongly in L 2 (il; M 3x3 ). Finally, as 

\ Xh v(h a - 1 G h )(G h ) T \ < h^Xh^ihi) \G h \ 2 < uj(h^) ff\G h \, 

also this last term converges to zero strongly in L 2 (Q,; M 3x3 ). Combining together 
(4.23) and the previous convergence properties, we obtain the claim (4.21). Notice 
that by (4.17) and (4.20) this implies E = CG e L 2 (fi;M 3x3 ). 

Step 4- Consequences of the Euler- Lagrange equations. We now begin to derive 
some preliminary information from the Euler-Lagrange equations (4.15). 

Let 4> £ C£(R 3 ;R 3 ) be such that (f>(x',x 3 ) = for every x £ Tx(-±, ±), and let 
us consider a test function of the form cj) h {x) := h<f>(x' , ^-). We notice that <j> h is 
an admissible test function, as (j) h & Cl(M. 3 ;R 3 ) and cj> h (x',hx3) = h(j){x' , X3) = 
for every x G rx(— | i). 
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Inserting (j> h in (4.15) leads to 

h f j2 RhEh ( Rh ) T ^- d A(y h )'^dx 

Jn i=l 

+ f R h E h (R h fe 3 -d 3( f > ((y h y, y f)dx= f tfge 3 ^{{y h )' V -f-)dx. 

As R h E h (R h ) T is bounded in L^ft; M 3x3 ) and V'<f> is a bounded function, the first 
integral on the left-hand side converges to zero as h — > 0. Since the right-hand side 
is clearly infinitesimal, we deduce 

Jim / R h E h (R h ) T e 3 ■ d 3 <f>((y h )' , dx = 0. (4.24) 
Jn " 

On the other hand, owing to (3.15), (4.5), (4.8), and to the continuity and bound- 
edness of d 3 cj), we have 

d 3 (j)((y h y, Hi) d S (f>(x',x 3 +v(x')) strongly in L 2 (S1;IR 3 ), if a = 3, (4.25) 

d^({y h )\ y) -» ^(a/.a*) strongly in L 2 (ft;R 3 ), if a > 3, (4.26) 

(the convergence is actually strong in L p (fi;R 3 ) for every p < oo). Therefore, 
splitting the integral in (4.24) as 

j^R h E h {R h ) T e 3 -d 3 J>((y h y y -j [ )dx 

= £ X hR h E h (R h fe 3 ■ d^{{y h )', ?±) dx 

+ £(1 - X h)R h E h (R h ) T e 3 ■ d 3 <t>{(y h )\ 
and using (4.20) and (4.21), we conclude that 

/ Ee 3 -d 3 (f>{x', x 3 + v(x'))dx = if a = 3, (4.27) 

[ Ee 3 -d 3 (j)dx = if a > 3, (4.28) 

Jn 

for every </> G C^(M 3 ;M 3 ) such that <j>(x',x 3 ) = for every x g Tx(-i, |). 

In the case a = 3, let Wk G C^M 2 ) be a sequence of functions such that the 
restriction of Wk to S converges to v strongly in L (S) and Wk{x') = for every 
x' G T. Then, given any <\> G Qj(M 3 ;K 3 ) satisfying </> = on rx(— \, i) wc can 
choose (j>k{x) := <j)(x' ,x 3 — Wk(x')) as test function in (4.27). Passing to the limit 
with respect to fc, we obtain that equation (4.28) holds true also for a = 3. 

From (4.28) it follows that Ee 3 = a.c. in ft. This property, together with the 
fact that E is symmetric, entails 

/ E n Ea \ 
E=\ E l2 E 22 (4.29) 
\ / 

for any a > 3. 
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Step 5. Zeroth moment of the Euler- Lagrange equations. Let E: S — > M 3x3 be 
the zeroth moment of the limit stress E, defined as 

E(x') := \ E(x)dx 3 (4.30) 



for every x' G S. In the following we derive the equation satisfied by E. 

We consider as test function in (4.15) a map independent of the variable x 3 . More 
precisely, let ip G (1R 2 ; IR 2 ) be such that ip(x') = for every x' G T. Choosing 
(j)(x) — (ip(x'),0) in (4.15), we have 

/ [R h E h (R h f]" : W((/)0 da = 0, (4.31) 

where [R h E h (R h ) T ]" denotes the 2x2 submatrix of R h E h (R h ) T , whose entries arc 
given by [R h E h {R h ) T \" j = R h E h {R h ) T ei ■ e h 1 < i,j < 2. 

As in the previous step, it is convenient to split the integral in (4.31) as 

f [R h E h (R h ) T ]" : Vip((y h )') dx = [ X h[R h E h {R h ) T ]" :V 4,{{y h )') dx 
Jn Jn 

+ f {l-Xh)[R h E h {R h ) T ]":V^{{y h )')dx. (4.32) 



By (4.5) and the continuity and boundcdncss of Vip, the sequence Vip((y h )') 
converges to V-0 strongly in L 2 (£l; M 2x2 ). Thus, by (4.21) we obtain 



lim / Xh[R h E h (R h ) T ]":Vii(y h Y)dx 



dx = I E": VV dx, 
n 



while, using the boundedness of Vifj and (4.20), we have that the last integral in 
(4.32) converges to 0, as h — > 0. Therefore, by (4.31) we conclude that 



/ E" : V'tp dx = 



for every ijj G C,J(K 2 ; M 2 ) such that = 0. In terms of the zeroth moment of the 
stress defined in (4.30), the previous equation yields 

E" -.V'tpdx' = (4.33) 

for every ip G C^(R 2 ;M 2 ) such that V|r = 0, and by approximation for every 
ip G H^S-^M. 2 ) with V|r = 0. 

Step 6. First moment of the Euler- Lagrange equations. We now derive the 
equation satisfied by the first moment of the stress, that is defined as 

t i 

E(x'):= J 2 x 3 E{x)dx 3 (4.34) 

for every x' G S. 

Let <p G Cj(M 2 ) be such that tp\r = and let us consider (j>(x) — (0, j^(p(x')) in 
(4.15). Since (4.5) and the continuity and boundedness of <p entail 

lim / gtp((y h )')dx = / gipdx= / gipdx', 
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we deduce that 

nl r 

g(pdx . (4.35) 



/ l J E[R h E h (R h ) T ]Md i <p((y h y)dx= f 

^ i— 1 S 



On the other hand, we can make partially explicit the limit on the left-hand side 
by applying the usual splitting O = B^ U (f2 \ B^) and considering the following 
decomposition: 

-R h E h (R h ) T = h a - 3 A h E h (R h ) T + h a - 3 E h (A h f + -E h , (4.36) 
h h 

where A h is the sequence introduced in (4.6). By the pointwise estimate (4.14) on 
E h and the fact that \A h \ < C/h a - 2 , we have 

\A h E h (R h ) T + E h (A h f\ < c( ^ (Id +^7 lG,i) + \A h \\G h \ 
so that, from the energy bound (3.12) we infer that 



f (1 - X h)\A h E h {R h ) T + E h (A h ) T \dx 
Jn 



I W{U ^ lGh) ^+\\G h h4(l-X h )A^ 



< Ch + C\\(l- X h)A h \\ 



L- 



As A h converges strongly in L 2 (S; M 3x3 ), we conclude that 

(1 - X h)(A h E h (R h ) T + E h (A h ) T ) ^0 strongly in i 1 (Q;M 3x3 ). (4.37) 

On the other hand, by (4.6) we have that A h — > A strongly in IP(S\ M 3x3 ) for any 
p < oo, while X hE h E weakly in L 2 (ft;M 3x3 ) by (4.21). Therefore, 

Xh (A h E h (R h ) T + E h (A h ) T ) -± AE + EA T weakly in L q {VL; M 3x3 ) (4.38) 

for any q > 2. Using (4.37), (4.38), and the fact that di(p((y h )') — > strongly in 
L p (il) for any p < oo, we deduce that 

f 2 

lim / y\A h E h {R h ) T + E h (A h ) T hidM(y h y)dx 



r 2 

/ y2[AE + EA T ] 3l di<pdx. (4.39) 
in 



In particular, from the definition (4.6) of A and the structure property (4.29) of E 
we have that [AE]3j = J2j=i EjidjV and [iM T ]3i = 0. Therefore, integrating over 
X3 we obtain 

f y^[AE + EA T ] 3i di^dx = f E" :(X7'v®V'ip)dx'. (4.40) 
Combining (4.35), (4.36), (4.39), and (4.40), we conclude that, if a = 3, 

2 

lim / J S^E% i d i ip{{y h )')dx = [ gtpdx' ~ f E" : (V v ® V <p) dx' (4.41) 
k->o in " ~i is is 

for every y> <E C^M 2 ) such that cp| r = 0. 
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If a > 3, by (4.39) we deduce that 

r 2 

lim / V h a - 3 [A h E h (R h ) T + E h (A h ) T ] 3i d lip {(y h )')dx = 0. (4.42) 
Combining this equation with (4.35) and (4.36), we conclude that, if a > 3, 

2 

, lim n / jJ2 E ^ d M(y h y)dx= [ g^dx' (4.43) 
Jn n l=1 Js 

for every <p G C£(R 2 ) such that <p|r = 0. 

We now want to identify the limit in (4.41) and (4.43) in terms of the first 
moment E. This will be done, as in [19], by first considering in the Eulcr-Lagrange 
equations (4.15) a suitable test function <j> with a linear behaviour in 23, and then 
passing to the limit with respect to h. Using the symmetry of E, this provides 
us with an identity (see (4.58) below) relating the first moment E with the limit 
in (4.41) (or (4.43)) and, by comparison with (4.41) (or (4.43), respectively) the 
limiting equation for E (see (4.62) and (4.63) below). 

An additional difficulty with respect to [19] is due to the fact that the simple 
choice 4>(x) = (2:377(2;'), 0) is not allowed in our framework, as this test function is 
not bounded in R 3 . For this reason we introduce a truncation function h . which 
coincides with the identity in an interval (~LOh,uJh), for a suitable u>h — * +00, and 
we consider a test function of the form (4.48) below. The rate of convergence of LOh 
has to be chosen in such a way to match two requirements. On one hand, we need 
to show that the limiting contribution due to the region where 9 h does not coincide 
with the identity is negligible. This can be done by means of the estimate (4.16), 
once we prove that the measure of the set Dh where \y^/h\ > u>h is sufficiently 
small. This is guaranteed if the rate of convergence of u>h is fast enough (see proof 
of (4.55) below). On the other hand, because of this choice, the L°°-norm of the 
test functions (j) h is not bounded, but blows up as u>h- Therefore, the convergence 
rate of Wh has to be carefully chosen to ensure that the integral on Q \ Bh remains 
irrelevant, as usual. This is possible owing to the choice of Bh and the estimate 
(4.19) (sec proof of (4.51) below). 

To be definite, let ujf, be a sequence of positive numbers such that 

huj h ^oo, hT-^Wh-yO, (4.44) 

where 7 is the exponent introduced in the definition of Bh- This is possible since 
7 < a — 2 (for instance, one can choose ojh '■= h~^ a t>' 2 ). Let 9 h E Cj(R) be a 
truncation function satisfying 

9 h (t) = t for\t\<w h , (4.45) 
\9 h (t)\ < \t\ for every t G R, (4.46) 

< 2. (4.47) 



L 3 



Let i] G C^(R 2 ; R 2 ) be such that 77(2/) = for every x' G T. We define 



H.S 



^ h (x):=(e h [^-)rj(x'),0). (4.48) 
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Owing to the assumptions on 9 h and rj, the (j> h are admissible test functions in 
(4.15); then inserting (j> h in (4.15) leads to 



+ f 1 j2[R h E'\RrUn l ((y h Y)(^(^))dx = o. (4.49) 

Jn i=1 

We now compute the limit of each term in (4.49) separately, starting with the first. 
We consider the usual splitting fl = B^ U (fi \ Bh) and we carefully analyse the 
contributions of the integral in the two subdomains. 
If a = 3, we have that 

lim / Xh 6 h (^-)[R h E h {R h ) T ]" :Vr)((y h y)dx= [ (E" + vE") : V'77 dx' . (4.50) 
h^o J Q \ h J J s 

Indeed, by (4.8) and (3.15) the sequence y^/h converges to x% + v a.e. in fl and is 
dominated by an L 2 function. From (4.45) and (4.46) it follows that the sequence 
(y^/h) converges to x$ + v a.e. in and is dominated by an L 2 function. Owing 
to the convergence (4.5) of y h and to the continuity and boundedness of V'77, we 
conclude that 

h (^)Vv((y h )') -» (x 3 + v)Vr)(x') strongly in L 2 (f!;R 2 ). 
Therefore, by (4.21) we deduce 

lim f Xh 6 h (^){R h E h (R h ) T }":V V ((y h y)dx= f (x 3 + v)E" : V V (x') dx. 
Jn v n ' Jn 

Integration with respect to x$ yields (4.50). 

As for the integral on Q\ Bh, by the estimate (4.47) on 9 h and (4.19) it can be 
bounded by 

| (1 - Xh)\0 h (^ ) [R h E h (R h ) T y : V'r,((y h y)\ dx 

< 2w h ||V^||i- I \E h \ < Ch^-^uh- (4.51) 
Jn\B h 

therefore, it is infinitesimal as h — > by the second property in (4.44). We conclude 
that, if a = 3, 

lim [ 6 h (^-)[R h E h {R h ) T }" ■.V'j 1 {{y h y)dx= [ {E" + vE") : Vrjdx'. (4.52) 
Jn v h ' Js 

Analogously, for a > 3, since y^/h converges to X3 strongly in L 2 (0), we deduce 
that 

lim f 9 h (^)[R h E h (R h ) T ]":V'r]{{y h y)dx= [ E":Vr]dx'. (4.53) 
Jn v h ) J s 
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In order to analyse the second integral in (4.49), it is convenient to split it as 
follows: 

oh . „,h s 



sE^*c*n.*«»*n(iU£))* 

2 — 1 



i=l 

+ /Jg [ ^ (flTla , i(( „7)(*(f)- 1 )^. (4.54, 

We claim that the second term on the right-hand side is infinitesimal, as h — > 0, 
that is, 

^ll± lR ^ { RrU^Yi^)- l) — 0. (455, 
JU i=l 

If a = 3, combining (4.49), (4.52), and (4.54), the claim implies that 

2 

lim f \y\R h E h {R h ) T ] iim {{y h )')dx = - f (E" + vE") : V^dx'. (4-56) 

On the other hand, we can identify the limit on the left-hand side of the equation 
above by applying the decomposition (4.36) and using the convergence properties 
(4.37) and (4.38); in this way, we deduce that, if a — 3, 

, lim n / \iyR h E h (R h ) T ]iMy h )')dx 

2 

= / E" : (r) ® Vv) dx' + lim f y V E% r H ((y h )') dx. (4.57) 
Js h ^°Jn h ~! 



Comparing (4.56) and (4.57), and applying equation (4.33) with tp = vi], we con- 
clude that, if a = 3, 

2 

i im n /lE ^3 *»((»*)') dx = — f E" : V V dx'. (4.58) 
If a > 3, combining (4.49), (4.53), (4.54), and the claim (4.55), we obtain 

2 

lim / \ Y\R h E h (R h ) T \ a m{{y h )')dx = - f E":V' V dx'; (4.59) 

moreover, by the decomposition (4.36) and the convergence properties (4.37) and 
(4.38) we have 

2 2 

, lim n / T t^[R h E h (R h ) T ]isVi((y h y)dx=]im [ ^Y,E^ m ((y h y)dx. (4.60) 
h ^°Jfi " ,- =1 h ^°Jn n i=l 

Comparing (4.59) and (4.60), we deduce that equation (4.58) is satisfied also for 
a > 3. 
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It remains to prove (4.55). To this aim we introduce the set Dh ■= {x G 17 : 
\y^(x)\/h > uJh}- Since the sequence y^/h is bounded in L 2 (fi) by (4.8) and (3.15), 
we have 



\D h \ I M d x<c^ x \D h \ l /\ 

JD h n 



which implies 



\D h \ < Cuj: 



(4.61) 



Since the derivative of 9 h is equal to 1 on (—uih,uJh) by (4.45), the integral in 
(4.55) reduces to 

/jD^^) T ]«"«^')(^(x)- 1 )^ 

By (4.16), (4.47), and (4.61), we have 

J Uh i=l 

C , 



"(' 



dt 

C, 



\vh- 



<Ch a - 2 + ^-\D h \ 1 / 2 <Ch a - 2 
h 



\E h \dx 



C 



huj h 



By (4.44) this proves the claim (4.55). 

Step 7. Limit equations. Let <p <G C 2 (M. 2 ) be such that <p(x') = 0, VV(s') = 
for every x' E T. Since E is symmetric, we can compare equation (4.58) (where we 
specify n = Vip) with (4.41), if a = 3, or (4.43), if a > 3. In this way we deduce 
that, if a = 3 



/ E" : (Vv ig> VV) dx' - [ E" :{V) 2 ipdx' = f g<pdx', 
Js Js Js 



while, if a > 3, 



E" :{V) 2 ipdx' = / gfdx' 



(4.62) 



(4.63) 



By approximation the two equations hold for every ip 6 H 2 (S) with <^|r = and 
Wlr = 0. 

In order to express the limiting equations (4.33), (4.62), and (4.63) in terms of 
the limit displacements, an explicit characterization of E" and E" is needed. Since 
E = CG and E is of the form (4.29), we have E" = C 2 G" (see [19, Proposition 3.2]). 
Therefore, by (4.10) and (4.11) we obtain, for a = 3, 

E" = £ 2 (symV'u + -V't; ® Vv), E" = —^£ 2 {V) 2 v. 

These identities, together with (4.33) and (4.62), provide us with the Euler-Lagrange 
equations (3.8)-(3.9). 
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By (4.10) and (4.12) we obtain, for a > 3, 

E" = £ 2 (symVy), E" = - J_£ 2 (V')V 

The first identity, together with (4.33) and the boundary condition u — on T, 
implies that u = 0, while the second identity, together with (4.53), provide us with 
the Euler-Lagrangc equation (3.11). This concludes the proof. □ 
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